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Abstract 

Emergence of hierarchies is investigated by Monte Carlo simulation in a timid 
society where all individuals are pacifist. The self-organiztion of hierarchies is 
shown to occur in two steps as the population is increased, i.e. there are three 
states, one egalitarian and two hierarchical states;the transition from the egali- 
tarian to the first hierarchical state is continuous and the transition from the first 
hierachical state to the second one is discontinuous. In the first hierarchical soci- 
ety, all individuals belong to either middle class or losers and no winners appear. 
In the second hierarchical society, many winners emerge and the population of 
the middle class is reduced. The hierarchy in the second hierarchical society is 
stronger than the hierachy in a no-preference society studied by Bonabeau et al 
[ Physica A217, 373 (1995)] 
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1 Introduction 

The emergence of hierarchies is a common phenomenon in societies and animal groups. 
In a pioneering work, Bonabeau et a/.£Q have shown that a hierarchical society can 
emerge spontaneously from an equal society by a simple algorithm of fighting between 
individuals diffusing on a square lattice. On the basis of results of Monte Carlo sim- 
ulation and an analysis by a mean field theory, they concluded that subcritical or 
supercritical bifurcations exist in the formation of the hierarchical structure as the 
density of individuals is varied. In their model, each individual is assumed to have 
some wealth or power which increases or decreases by winninng or losing in a fight. 
The essential processes of the model are diffusion, fighting and spontaneous relaxation 
of the wealth. Various societies can be modelled by specifying each process and the 
emergence of the hierarchy depends strongly on the specifications. [21 El 

In this paper, we investigate a variation of the model introduced by Bonabeau et 
al. P, where the diffusion algorithm is modified to include the effect of the trend of 
society. Namely, we study the emergence of hierarchies in a timid society, in which 
an individual always tries to avoid fighting and to fight with the weakest among the 
neighbors if he/she cannot avoid fighting. By Monte Carlo simulation, we show that 
the emergence of the hierarchy is retarded in the timid society compared to the no- 
preference society investigated by Bonabeau et al. and that the transition to the 
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hierarchical state occurs in two successive transitions of a continuous and a discontin- 
uous ones. Consequently, there exist three different states in the society, one equal 
and two hierarchical states. In the first hierarchical states, we see no winners but 
losers and people in the middle class. In the second hierarchical states, many winners 
emerge from the middle class. We also show that the distribution of wealth in the sec- 
ond hierarchical state is wider compared to the hierarchical state of the no-preference 
society. 

In Sec. 2, our model is explained in detail. Results of Monte Carlo simulation are 
presented in Sec. 3. In Sec. 4 the characteristics of the hierarchical states is analyzed 
in detail. Section 5 is devoted to discussion. 



2 A timid society 



We consider N individuals diffusing on an L x L square lattice, where every lattice 
site accomodates at most one individual. An individual is to move to one of nearest 
neighbor sites according to the following protocol. When individual i tries to move to a 
site occupied by j, i and j fight each other. If i wins, i and j exchange their positions, 
and if i loses, they keep their original positions. We associate each individual a quantity 
which we call power or wealth. The power increases by unity for every victory and 
decreases by unity for every loss. The probability Qij that % wins the fight against j is 
determined by the difference of their powers Fj and Fj as 

QlJ = l + exp^-F,)]' (1) 

where rj is introduced as a controlling parameter. When rj = oo, the stronger one 
always wins the fight and when 77 = 0, the winning probability of both ones are equal. 
We also assume that the power of individuals relaxes to zero when they do not fight, 
namely power Fi(t + 1) at time t + 1 is given by Fi(t) throughp^ 

Fi(t + 1) = Fi(t) - /itanh[Fi(*)]. (2) 

Here, the unit of time is defined by one Monte Carlo step during which every individual 
is accessed once for move and \i represents an additional controlling parameter. This 
relaxation rule indicates that people lose their wealth of a constant amount when their 
power is large, and when their power is small, they lose it at a constant fraction, 
namely they behave rather miserly. It also indicates that the negative wealth (debt) 
can relax to zero in the similar manner. Note that this rerlaxation rule is critical to 
the emergence of hierarchical society. |2j 

We characterize the timid society by the preference of individuals in diffusion. In 
the timid society, every individual favors not to fight and thus it moves to a vacant 
site if it exists. If no vacant sites exist in the nearest neighbors, then it moves to a site 
occupied by an individual whose power is the smallest among the neighbors. When 
more than two neighbors have the equal power, then an opponent is chosen randomly 
from them. 

We characterize the static status of the society by an order parameter a which is 
defined by[HI2] 
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Here, N is the number of individuals, and Di and Si are the number of fights won and 
lost, respectively, by individual i. Note that o = corresponds to an egalitarian status 

and cr = l/\/l2 ~ 0.2887 when the chance for victory jjrqj c;. is distributed uniformly 
in [0, 1]. After sufficiently long Monte Carlo simulation, variation of a is stabilized and 
one can use it as an order parametrer. 

We also monitor the population profile by focusing on the winning probability. We 
classify individuals into three groups by the number of fights which an individual won; 
winners are individuals who won more than 2/3 of fights and losers are individuals who 
won less than 1/3 of fights. Individuals between these groups are called middle class. 

3 Monte Carlo simulation 

Monte Carlo simulation was performed for N = 3500 individuals on the square lattice 
with periodic boundary conditions from L = 60 to L = 180. We obtained the order 
parameter a 2 and other quantities for 10 6 Monte Carlo steps. 
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Figure 1: Order parameter a 2 as a function of p = N/L 2 with p = 0.1, for four different 
values of rj: rj — 50(D), 0.5(x), 0.05(0)) 0.005(A). Error bars are much smaller than 
the size of symbols. 

Figure 1 shows the dependence of the order parameter a 2 on the density p = N/ L 2 
for several values of rj, where p is fixed to p = 0.1. We can see two clear transitions; 
one at a lower critical density pci and the other at a higher critical density pc2- 
The transition at pci is continuous and the transition at pc2 is discontinuous. The 
dependence of the critical densities pci and pc2 on parameter rj is shown in Fig. 2. 

We can identify three states for a given value of rj; an egalitarian state for p<pci, 
a hierarchical society of type I for pci<P<Pc2 and a hierachical society of type II for 
Pc2<P<X- I n the egalitarian society, winners and losers lose their memory of previous 
fight before they engage in the next fight, and thus they changes their status in time. 
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Figure 2: The dependence of the critical densities pci (the circles) and pc2 (the crosses) 
on parameter rj, where p — 0.1. The curves are the guide for eyes. 

In the hierarchical state, a winner keeps winning and a loser keeps losing. We discuss 
the difference between type I and type II hierarchical society in the next section. 

The results strongly depend on p. We show the phase boundary on the p-p plane 
for 7] = 0.05 in Fig. 3. 

4 Two hierarchical societies 

In order to investigate the structure of the hierarchical states, we analyze profile of 
population. The dependence of the population of each class is plotted against the 
density in Fig. 4. Rapid changes of the populations signify emergence of different state 
of the hierarchical societies. In the egalitarian state p<pci, all individuals belong to 
the middle class as expected. In the hierachical society I pci<P<Pc2, some individuals 
become losers whose number increases as the density is increased, but no winners are 
seen. In the hierachical state II p > pc2, many winners appear and the population in 
the middle class is reduced significantly. 

Figures 5 (a), (b) and (c) show the spacial distribution of individuals after 10 6 
Monte Carlo steps in the egalitarian, the first hierarchical and the second hierarchical 
states, respectively. No specific spatial inhomogeneity is observed in the timid society. 

In order to see details of the hierarchical structure, we plot the population as a 
function of the density and the winning probablity in Fig. 6. From this plot, we 
conclude that (1) in the hierarichical state I, people in the middle class with slightly 
higher winning probability increase, but no winners are seen and (2) in the hierarchical 
state II, the most of winners have very high winning probability, while people in the 
losers and the middle class are distributed in a wide region of the winning probability, 
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Figure 3: The dependence of the critical densities pci (the circles) and pc2 (the crosses) 
on parameter \i for rj = 0.05. The curves are the guide for eyes. 

5 Discussion 

We have investigated the emergence of self-organized hierarchies in the timid society. 
Our results show that the emergence of the hierarchical state in the timid society is 
retarded compared to the no-preference society. This delay is natural since individuals 
in the timid society tend to avoid fighting and thus the wealth is distriruted more or less 
evenly among individuals when the population is low. Furthermore, the emergence of 
the hierarchical society in the timid society occurs in two steps, and the first transition 
is continuous and the second one is discontinuous. The strength of the hierarchy in 
the high density region is stronger in the timid society compared to the no-preference 
society. For the same choice of t] = 0.05 and \i — 0.1, a 2 for the former case is twice as 
large as the latter pQ. 

To understand these behaviors, we first remind the fact that the the hierarchical 
society emerges when the power cannot relax before the subsequet fight. In the timid 
society, an idividual can avoid fighting when the density is low, and thus the ealitarian 
state is favored for low densities. In the timid society, weaker individuals have more 
chance to be challenged and thus to lose their power, and stronger ones has less chance 
to fight and their power stay near zero. This situation corresponds to the hierarchical 
state I. When the density is increased above the upper critical density, all individuals 
have more chance to fight and thus stronger individuals become much stronger. 
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Figure 4: Dependence of the population in each class on the density when /i = 0.1 and 

T] = 0.1. 
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Figure 5: Structure in the quilibrium state, (a) the eagalitarian state p = 0.3, (b) the 
hierarchical socity I p = 0.5 and (c) the hierarchical socity II p = 0.7. the circles are 
the winner, the traiangles are individuals in the middle class and the crosses are the 
loser. 
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